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Based on scanning tunneling microscopy and first-principles theoretical studies, we characterize 
the precise atomic structure of a topological soliton in In chains grown on Si(lll) surfaces. Variable- 
temperature measurements of the soliton population allow us to determine the soliton formation 
energy to be ~60 meV, smaller than one half of the band gap of ~200 meV. Once created, these 
solitons have very low mobility, even though the activation energy is only about 20 meV; the sluggish 
nature is attributed to the exceptionally low attempt frequency for soliton migration. We further 
demonstrate local electric field enhanced soliton dynamics. 

PACS numbers: 73.20.Mf, 68.35.Md, 68.37.Ef 



One-dimensional (ID) materials have attracted inten- 
sive attention due to their richness in physics and poten- 
tial applications in nanoelectronics. The significantly en- 
hanced interactions between charge, spin, and lattice in 
such ID electronic systems could lead to exotic ground 
states, such as Luttinger liquid [lj], and charge density 
waves (CDWs) 0. As prototype model systems, metals 
on semiconductor surfaces at low coverages could be self- 
organized to form atomic chains Q, which may exhibit 
ID CDW states at low temperatures Q. The underly- 
ing physical properties of such collective electronic states 
in ID systems, in turn, could be elucidated with atomic 
precision, for example by scanning tunneling microscopy 
(STM). 

Besides the presence of exotic ground states, poten- 
tially even more intriguing are the elementary excitations 
of the ID CDWs. Examples include the excitations of 
phase and amplitude (known as phason and amplitudon 
[5(), and the nonlinear topological excitation or soliton, 
which has been well studied in ID conjugated polymers 
such as polyacetylene @. A soliton can be regarded as 
a local phase-flip boundary separating two energetically 
degenerate ID CDW phases. The lattice displacements 
from the undistorted x 1 structure occur in opposite di- 
rections across a single soliton, thereby giving rise to 
a geometrical feature that is more complex than a sin- 
gle phase boundary. The formation energy of a soliton 
should also be less than one half of the band gap, oth- 
erwise electron-hole pairs would be more readily excited 
upon increasing the temperature. Solitons may possess 
spin-charge inversion properties, and act as the effective 
carriers that account for the high conductivity in con- 
ducting polymers @. However, earlier studies of solitons 
were based on ensemble average techniques Q , and only 



very recently had some preliminary qualitative reports 
been made on the existence and dynamics of solitons us- 
ing STM 0-i]. To date, a comprehensive quantitative 
study of topological solitary excitations at the atomic 
level remains a challenge. 

In this Letter, we present a comprehensive study of 
the structure, energetics, and dynamics of topological 
solitary excitations at the atomic scale using a combina- 
tion of STM and first-principles techniques. The precise 
atomic structure of the solitons in In chains on Si(lll) 
surface is determined reliably. The formation energy and 
diffusion barrier of the solitons are estimated by examin- 
ing the soliton population and thermal dynamics at dif- 
ferent temperatures. We further demonstrate the capa- 
bility of manipulating the soliton dynamics by applying 
an electric field between the STM tip and the solitons. 
Together, our findings represent the first insight into 
atomic-scale characteristics of topological excitations in 
such ID systems. 

The experiments were performed in an ultra-high vac- 
uum system, equipped with an Omicron low temper- 
ature (LT) STM with a base pressure below lxlCT 10 
mbar. The Si(lll) substrates (n-type, ^0.001 Vt-cm) 
were cleaned using the standard "flashing" recipe until 
the 7x7 reconstruction was established at the surface. 
Then, the Si(lll)-(4x 1)-In atomic chain structure was 
prepared by evaporating about one monolayer In onto a 
clean Si(lll) surface at 700 K, followed by postannealing 
at the same temperature for 30 min. 

As shown in Fig. 1(a), the STM image of the In chains 
on Si(lll) at room temperature (RT) clearly shows the 
4x1 pattern. As the temperature is reduced, the system 
undergoes a reversible phase transition from a metallic 
4x1 phase to an insulating 8x2 phase at a transition 
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FIG. 1. (color online) Typical STM images of In chains on 
Si(lll) surface at RT (a) and 78 K (b) at sample bias (V s ) 
of -1.1 V and 1.5 V, respectively. The 4x1 and 8x2 unit cells 
are drawn in (a) and (b), respectively. The inset of (b) shows 
the two 4x2 phases with alternating orientations composing 
the 8x2 structure, (c) dl/dV spectroscopes taken at the 8x2 
region, (d) Experimental STM images of the same soliton at 
V s of 0.6 V (top), 1.5 V (middle), and -0.6 V (bottom), (e) 
Simulated STM images of a soliton at V s of 0.6 V, 1.5 V, 
and -0.6 V from top to bottom, respectively. The adopted 
optimized soliton structure model is shown at the bottom. 



temperature around 130 K 10]. This is evidenced by 
the 8x2 pattern in the STM image at 78 K (Fig. 1(b)). 
However, the structural formation of the 8x2 phase is 
still controversial: One possible scheme is the pairing of 
the outermost In atoms, resulting in the formation of In 
the other includes additional shear distor- 
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tion between two neighboring dimerized chains, leading 
to the formation of In hexagons [12( . From previous den- 
sity functional theory (DFT) calculations, it is found that 
the former structure is still metallic fl3| | . while the latter 
is insulating 12j. Here we acquired dl/dV spectroscopy 
on the 8x2 surface by adopting a locking- in technique, 
and the typical result is shown in Fig. 1(c). A fully- 
opened energy gap of about 0.2 eV is clearly observed, 
which supports the hexagon structural model. 

The 8x2 structure is composed of two 4x2 chains with 
alternating orientations denoted as A and B phases in the 
inset of Fig. 1(b). The two 4x2 phases are degenerate in 
energy with a gliding reflection symmetry It is therefore 
possible that the A and B phases coexist in the same 
4x2 chain, and produce a phase boundary, namely, a 
topological solitary excitation between them. 

In the LT STM images [see Fig. 1(b)], we find 
some bright features with well-defined uniform structure, 
which are absent at RT. More importantly, the orienta- 
tions in the In chains are mirror-symmetric around the 
bright features, i.e., the bright features act as phase- flip 
boundaries separating the energy-degenerate A and B 
phases. Zoom-in STM images at different V s are shown 



in Fig. 1(d)- The 4x2 structure shows Vs-dependent 
patterns [l4| , which are always mirror-symmetric around 
the phase boundaries. These phase boundaries are lo- 
calized and preserve their structure even after diffusion, 
thus are likely to be originated from the formation of 
solitons. 

In order to determine the precise soliton structure, we 
perform first-principles DFT calculations [15| within the 
generalized-gradient approximation [la ]. All atoms are 
described by norm-conserving pscudopotcntials 17| . On 



the basis of the 4x2 hexagon model [l2[ , we simulate the 
soilton structure by using a periodic slab geometry, where 
each slab contains six Si atomic layers (not including the 
Si surface chain) and the bottom Si layer is passivated by 
one H atom per Si atom. The vacuum spacing between 
these slabs is 10 A. The electronic wave functions are 
expanded in a plane- wave basis set with a cutoff energy 
of 15 Ry. The k-space integrations are done with four 
points in the surface Brillouin zone of the 4x17 unit cell. 
Here we consider the phase boundary by employing the 
4x17 supercell where two solitons (separated by three In 
hexagons) are included because of the use of a periodic 
slab geometry. Based on our experimental evidence, we 
consider mirror symmetric In hexagons with respect to 
a plane of the phase boundary. The optimized soliton 
structure with a heart-shaped phase boundary is shown 
in Fig. 1(e). 

For comparison with the STM observations, we simu- 
late the constant-current STM images using the Tersoff- 
Hamann approximation [l8j , all at the charge density of 
p = 10~ 6 electrons/bohr 3 . The results are displayed in 
Fig. 1(e). We find that for all the bias voltages con- 
sidered, the orientations of the A and B phases as well 
as the separation of the bright spots agree very well be- 
tween the theoretical simulations and experimental data. 
As shown in Fig. 1(e), the bright spots composing A 
and B orientations can be identified in terms of the two 
outer In atoms. It is interesting to note that, when the 
bias voltage changes from 0.6 V to 1.5 V, the orienta- 
tions of A and B phases are swapped with each other in 
Fig. 1(d), which is caused by the slight lateral shifts of 
the bright spots between the 0.6 and 1.5 V images. The 
simulated images with V s = 0.6 and 1.5 V also repro- 
duce well the features at the phase boundary as shown 
in the experimental ones. However for V s ~ -0.6 V, the 
simulated image at the phase boundary differs somewhat 
from the experimental image. This discrepancy is likely 
to be caused by the apparently different features in the 
simulated filled-state image, where weak protrusions orig- 
inating from inner In atoms exist between adjacent paired 
bright spots, whereas the corresponding experimental im- 
age is unable to resolve such weak protrusions. Namely, 
the discrepancy is primarily due to the resolution limita- 
tion of the STM tip in resolving dark features [lj| . 

Temperature-dependent experiments are performed to 
explore the energetics of the solitons, with temperatures 
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FIG. 2. (color online) (a) Temperature-dependent soliton 
density n during warming-up. The large filled circle denotes 
n at 78 K after cooling-down, (b) Temperature-dependent 
hopping rate v s . The inset of (b) shows the time-dependent 
probability p that solitons do not move during scanning at 86 
K. 



ranging from 78 K to 130 K. At temperature above ~110 
K, we observe the coexistence of the 4x2 (8x2) and 4x1 
phases [20] , and focus on the soliton behavior in the 4x2 
(8x2) phase. In determining the soliton formation en- 
ergy, we measured the soliton density n at various tem- 
peratures using the following procedure: At a given tem- 
perature, several well-defined 4x2 (8x2) regions were 
scanned, and the total number of isolated solitons was 
accumulated. To minimize the fluctuations in the soli- 
ton number, we ensured that the total soliton number 
counted at each temperature be about 500. The soliton 
density n was then obtained by taking the total soliton 
number divided by the total area scanned. We checked 
the thermal reversibility by measuring the soliton popu- 
lation first by warming up and then by cooling down the 
systems. For the cool-down process, we only counted the 
soliton density at 78 K after cooling-down, and its value 
is close to that before warming-up as shown in Fig. 2(a). 
This cross check convincingly demonstrated that the soli- 
ton density does not depend on the thermal history, and 
the solitons (aside from that associated with uq discussed 
below) are indeed thermally excited. 

Thus, the formation energy E c for a soliton can be 
evaluated by fitting the data using the formula n — 
n +niexp(-i? c / k B T). We obtain E c of 61 ± 14 meV, less 
than one half of the band gap (200 meV). It means that 
the solitons, rather than the electron-hole pairs, are fa- 
vored to be thermally excited [f| , consistent with the ob- 
servations that solitons are stabilized in the In/Si(lll) 
system. This further supports the soliton formation. 
Here the residual soliton density no is estimated to be 
about 0.06 nm -2 . As defects are unavoidably present 
on the surface, some of the thermally excited solitons 
close to the defects could be stabilized to become immo- 
bile, and will not annihilate during cooling either, thereby 
leading to a non-zero uq. This observation has been fur- 
ther verified by going to 5 K in our STM observations: at 
this exceptionally low temperature, solitons can still be 
observed with a density close to no. Using DFT calcula- 
tions, we estimate the formation energy and the energy 



gap to be about 110 and 200 meV, respectively, close to 
the experimental values. In contrast to the experimental 
data, the calculated formation energy is slightly larger 
than half of the energy gap. This may be due to the 
well-known bandgap underestimation of the DFT calcu- 
lations. 

The importance of soliton migration lies in the aspect 
that solitons act as the effective carriers for high conduc- 
tivity in conducting polymers. It is likewise intriguing to 
examine the soliton dynamics in the surface-based atomic 
chains scale as well, as such chains may serve as elemental 
building blocks in future nanoelectronics. From the time- 
dependent probability p(t) that solitons remain immobile 
during successive STM scanning, we can determine the 
hopping rate v s at a given temperature by fitting p(t) 
with p(t) = po + (l-po)exTp(-i> s t) . As an example, the fit- 
ting at 86 K is shown in the inset of Fig. 2(b). The non- 
zero po is associated with no discussed earlier, namely, 
some solitons close to pre-existing defects are trapped to 
be immobile. By fitting the temperature-dependent v s 
with the Arrhenius equation v s = Aexp(-E a /kBT) [2l|, 
the diffusion barrier E a and the effective attempt fre- 
quency prefactor A are estimated to be about 16 meV 
and 6xl0 -3 Hz, respectively. 

Here we note that for typical surface diffusion, there 
exists an effective activation energy measuring the diffu- 
sion barrier, no matter what kind of diffusion mechanism 
involved. Usually the diffusion still follows the Arrhe- 
nius law, and the examples even include long jumps [22} 
and defect-mediated diffusion (23|. It is reasonable to fit 
the Arrhenius law to extract the activation energy in our 
case, and the excellent fitting displayed in Fig. 2(b) fur- 
ther confirms that an overall effective activation barrier 
against solition migration is well defined. 

We note that anomalously low attempt frequency of 
the order of 10 3 Hz has been reported for Al atom dif- 
fusion on inhomogeneous Au(lll) surface, which is as- 
sociated with the low activation energy (30 meV) [24$ . 
Generally observed trend is that the attempt frequency 
down shifts tremendously with shrinking diffusion barrier 
when the diffusion barrier is smaller than 100 meV [HI, a 
phenomenon commonly known as the kinetic compensa- 
tion effect [23. More importantly, unlike single atom or 
molecule, here the soliton is a collective topological exci- 
tation. Its diffusion must involve the rearrangements of 
many local bonds between the constituent atoms defin- 
ing the soliton. Therefore, the diffusion of a soliton is 
in principle a collective behavior, and thus may possess 
an extremely low attempt frequency. We expect that 
the present quantitative study on the soliton mobility, 
with exceptionally low activation energy and attempt fre- 
quency, will stimulate future efforts on trying to uncover 
the precise atomistic mechanism for solition migration. 

We further explore the feasibility to enhance the soli- 
ton migration by applying a voltage pulse between the 
STM tip and the solitons at 78 K [27J. The experimen- 
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FIG. 3. (color online) The displacement probability per volt- 
age pulse P dependent on the pulse voltage. The inset is the 
I-t curve during a pulse with V s = 1.6 V. 



tal procedure is as follows: The STM tip is scanned over 
the surface with V s = 0.8 V and / = 0.5 nA. Then it is 
moved to just over a soliton and a sample bias pulse with 
a period of 3.6 s is applied after disabling the feedback. 
The inset of Fig. 3 displays a recorded I-t curve when a 
pulse with V s — 1.6 V is applied. A steep decrease of the 
tunneling current is observed, indicating that the soliton 
moves away from its original position during the pulse. 

Here we found that the displacement probability per 
pulse P has a linear dependence on the bias voltage with 
about 80-100 events counted at each bias. Furthermore, 
there exists a clear polarity asymmetry: The displace- 
ment probability is much higher for positive bias than 
that for negative bias. For example, P is almost 0.8 for 
a pulse with V s = 3.0 V, while only 0.2 for a pulse with 
V s = -3.0 V. This strong polarity asymmetry is possibly 
due to the enhanced Coulomb repulsion between the ex- 
tra electrons and the n-type substrate when the electrons 
are injected into a soliton. It is noted that similar mani- 
festation of the substrate doping effect on the electronic 
properties of surface layers has been previously reported 
[28} . These experimental findings clearly demonstrated 
the capability of enhancing the soliton dynamics at the 
atomic scale. 

We may also acquire the effective mass of a soliton by 
adopting the standard procedure 0. The soliton length 
21 (in unit of a) is estimated to be roughly 5, from the 
full-width (along the chain direction) at half-maximum 
of the soliton feature in the STM image with V s = 1.5 
V (after filtering the x2 period modulation). The ef- 
fective mass of the soliton can be obtained from M s = 
(4/3Z)(urj/ a) 2 M, where uq is the lattice distortion of the 
x2 phase relative to the xl phase, and M is the mass 
of the In4 group. With uq = 0.35 A (estimated from 
the optimized structure) and I — 2.5, M s is assessed to 
be about 3700 m e (m e denotes the electron rest mass), 
about twice of the mass of a hydrogen atom. 

In polyacetylene, solitons have long length (21 — 14), 
low diffusion barrier (2 meV), and light effective mass 
(6 m e ), acting as excellent carriers [6j. In contrast, the 
solitons in the In/Si(lll) chains possess shorter length, 
higher effective mass, higher diffusion barrier. Together 
with the extremely low effective attempt frequency, they 



account for the sluggish nature of the solitons in the 
In/Si(lll) chains. 

In summary, we have characterized the precise atomic 
structure of a topological soliton in In chains grown 
on Si(lll) surfaces, by combining STM measurements 
and theoretical simulations. The soliton formation en- 
ergy is estimated to be ~60 meV from the temperature- 
dependent soliton population, which is smaller than one 
half of the band gap. The solitons possess very low mo- 
bility, even though the activation energy is only about 
16 meV; the sluggish nature is attributed to the excep- 
tionally low attempt frequency for soliton migration. Lo- 
cal electric field enhanced soliton dynamics has been fur- 
ther demonstrated. Our work gains the first insight into 
characteristics of topological solitary excitations at the 
atomic scale in ID electronic systems 
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